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OO ■ Abstract The quasifree reaction np — * pp7[~ was studied in a kinematically complete 

o ■ 

I experiment by bombarding a liquid hydrogen target with a deuteron beam of momentum 



1.85 GeV/c and analyzing the data along the lines of the spectator model. In addition to the 
three charged ejectiles the spectator proton was also detected in the large-acceptance time- 



H ' of-flight spectrometer COSY-TOF. It was identified by its momentum and flight direction 
thus yielding access to the Fermi motion of the bound neutron and to the effective neutron 
4-momentum vector P„ which differed from event to event. A range of almost 90 MeV excess 
energy above threshold was covered. Energy dependent angular distributions, invariant mass 
spectra as well as fully covered Dalitz plots were deduced. Sizeable pp FSI effects were found 
as were contributions of p and d partial waves. The behavior of the elementary cross section 
(Joi close to threshold is discussed in view of new cross section data. In comparison with 
existing literature data the results provide a sensitive test of the spectator model. 
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1 Introduction 



The low-energy regime in strong interaction physics is known to be dominated by pion 
producing reactions. Pioneering experiments were first carried out with bubble chambers 
which somewhat later were superseded by more refined experiments with electronic detectors 
at such hadron facilities as TRIUMF, LAMPF, PSI, and SATURNE and then, starting in the 
late 1980s, at the cooler rings lUCF, CELSIUS, and COSY. A vast amount of data on meson 
production in general and on pion production in particular became available, as can be seen 
in the recent review articles from the experimental [1] and the theoretical side [2]. Assuming 
isospin invariance, all single pion production reactions in A^A^ collisions with three body final 
states can be decomposed into a sum of at most two out of three elementary cross sections 
c/i,// [S] with Jj and If denoting the isospin of the A^A^ system in the initial and final state, 
respectively. Starting from the pp-entrance channel two reactions are possible, pp — > ppir^ 
and pp puTT^, yielding, respectively, an and an + aio. The isoscalar cross section aoi 
can only be determined in neutron-proton collisions as, e.g., anp^ppn- = + o"oi) and 

hence 

C"01 2 ■ (Tjip—^ppTY- ^pp^ppTT^ • (1) 

Close to threshold only a few partial waves have to be considered. In a simple, semiclassical 
picture the maximum / = Imax is given by Imax — R-Q- Here R, which gives the distance from 
the center of the A^A^ collision where a pion of momentum q is created, is of order h/ {m-i^c), 
i.e., the compton wavelength of the pion, and hence Imax — Q/^n — V- This dimensionless 
parameter rj is often found in the literature for specifying excess above threshold. A more 
direct measure is the excess energy Q = y/s — ^/so with ^/s{^/so) denoting the total CM 
energy (the threshold energy), respectively. Since for excess energies less than 100 MeV r] 
is well below 1.5, only partial waves of type Ss,Sp,Ps and Pp as, e.g., ^Pq ^^SqSq have 
to be considered, whereas the role of / = 2 contributions {Sd and Ds) can be regarded as 
negligible. Here the Rosenfeld notation Lplg has been used [3] with Lp being the orbital 
angular momentum of the A^A^ pair, Iq that of the pion with respect to this pair. 
In the past the pp entrance channel attracted most of the attention. Only recently the focus 
was shifted to a larger degree towards the study of proton-neutron collisions. In this case one 
might either use a neutron beam impinging on a hydrogen target or use a deuterium target 
as a substitute for a neutron and use the spectator model to determine the observables of 
the quasifree reaction. The former approach was extensively used by the Freiburg group 
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working at PSI [4-7] and to some extent by groups at LAMPF [8] and TRIUMF |9]. The 
first dedicated experiments on a deuterium target were also performed at TRIUMF [10,11] 
wfiere in order to bypass tlie difficulties of low-energy spectator detection, the final-state 
protons were restricted to the ^Sq, the "diproton" state, by selecting only two-proton events 
with a small relative momentum. Alternatively one can use a deuteron beam hitting a 
hydrogen target. The advantage of having fast spectator protons flying in the forward 
direction is counterbalanced by the fact that appropriate deuteron beams beyond 4GeV/c 
are not available any more and thus one is restricted to reactions where only light mesons 
are produced. The feasibility of this approach was investigated in our recent paper [12] on 
the dp —* pppn~ reaction. 

It is the aim of the present paper on the quasifree np pp'K~ reaction to present angular 
distributions, Dalitz plots and invariant mass distributions for the three reaction products, 
from which one might deduce clues as to the participating partial waves. We will present our 
data for angles as defined in Fig. 1. In the 3-particle center-of-mass (CM) system all three 
momentum vectors lie in a plane, their sum adds up to zero. Observables when given in the 
CM system will be marked with an asterisk. Denoting the beam direction by A^, the pion 
momentum by q* and the two proton momenta by p*, we will present angular distributions 
in the two angles cos 6* = cos (g*, A^) and cos 6p = cos (P*, A^) where P*, the relative proton 
momentum, is taken as the difference vector P2 By construction the latter is symmetric 
around cos 9p = 0. Strong interference effects between outgoing Ss and Sp on one hand and 
Ps and Pp waves on the other will produce highly asymmetric distributions in cos 6* even 
at the lowest Q(<20MeV). 
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2 Experimental Procedure 



2.1 Principle of measurement and spectator tagging 

The reaction under study is np ppn^ . Due to the lack of a neutron beam, deuterons 
were used instead and the data were analyzed along the lines of the spectator model. This 
method has been described in detail in [12], hence we will give only some short remarks here. 
The basic idea of the model is that 1) the proton in the deuteron can be regarded as an 
unaffected spectator staying on-shell throughout the reaction and 2) the matrix element for 
quasi-free pion production from a bound neutron is identical to that for free pion production 
from an unbound neutron. Crucial to the method is the task to detect and identify the 
spectator proton ps, since the information gathered from this particle gives a direct measure 
of the Fermi momentum carried by the off-shell neutron within the deuteron at the time of 
the np reaction. The Fermi momentum distribution as calculated from any of the existing 
A^A^ potentials has a maximum near 40 MeV/c and a tail extending towards several hundred 
MeV/ c, hence a wide range in excess energy Q can be covered with a monoenergetic deuteron 
beam. The main result of our recent study was that the two assumptions quoted above can 
be regarded as being fulfilled for Fermi momenta below 150MeV/c. 

The experiment was carried out with the time-of-flight spectrometer COSY-TOF set up 
on an external beamline of the proton synchrotron COSY [13] at the Forschungszentrum 
Jiilich. A deuteron beam of momentum pd = 1.85 GeV/c was focussed onto a liquid hydrogen 
target, charged particles emerging from the reaction zone were detected in a multi-layer 
scintillator hodoscope with its main components Quirl, Ring and Barrel. Details of the 
various subdetectors, their performance as well as the different steps necessary for calibrating 
the whole system have been described in a series of papers, see [12] and references therein. 
Here only a short overview will be given. By measuring each particle's flight time and 
direction, their velocity vectors given as v = {f3, 6, 0) could be determined with a time-of- 
flight resolution of better than 300 ps (cr) and an angular track resolution of better than 
0.3° (a). The momentum 4- vectors P of all detected particles were then obtained from 
the measured observables by applying additional mass hypotheses. Carrying out various 
tests as, e.g., momentum conservation, missing mass and invariant mass analyses as well 
as comparisons with results obtained from our Monte Carlo simulations helped to find the 
correct assignment for each event with a high degree of probability as quantified below. 

In the reaction dp ppir'ps four charged particles are emitted which in most cases all are 
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detected in the time-of-flight spectrometer. Thus the main trigger condition was such that 
a total of four hits was required in any of the stop scintillator hodoscopes Quirl, Ring and 
Barrel and at least one hit in the twelve-fold segmented start scintillator. Due to the fact 
that pions can also be emitted into the backward region where no detector was installed we 
set up a second trigger condition with only three required hits at a reduction factor of 10. 
Since for these events the unobserved pion can be reconstructed through a missing mass 
analysis, the full kinematically allowed phase space was covered. 

With a beam momentum below the threshold for 27r-production, any 4-hit event apart from 
accidentals could only result from the reaction under study. As the first step in our analysis 
we checked on the four possible hypotheses, i.e., the pion being particle 1, 2, 3, or 4 and 
calculated for each case the sums of longitudinal and transversal momentum components 
YIpl and YIpt- As the correct assignment we took the one where these values were closest to 
J2PL = Pd and J2pt = 0- As the spectator proton we then chose the one which was detected 
close to the beam axis with a momentum near Pd/2. The spread in Fermi momentum caused 
the momentum of the spectator to vary considerably, higher Q values correspond to lower 
momenta and vice versa. This is illustrated in Fig. 2 where for two narrow ranges in Q 
([18.0-34.0] and [61.0-74.5] Me V) the momentum distribution of the spectator proton and 
the summed distribution of both reaction protons is plotted. The spectator distribution 
given by the solid histogram at a mean < Q > = 26 MeV sticks out as a sharp line well 
separated from the much broader momentum distribution of the other two protons, whereas 
at < Q > = 68 MeV the spectator line is still rather narrow, but starts to overlap with the 
one of the reaction protons. Since a unique identification of the spectator is essential for 
the analysis we found it necessary to also limit the range in Q due to this effect and only 
considered events where the excess energy was below 90 MeV which roughly coincides with 
our proposed limit for the Fermi momentum p2]. In our finally accepted data sample of 
2.2-10^ events we obtained a longitudinal momentum distribution Ypl which had its center 
at pd with a width of 39MeV/c (a). In case of the transversal distribution the spread was 
even smaller, namely 13MeV/c. Alternatively we also used the missing mass method for 
identifying the various ejectiles, see ref.[12], and found full agreement. 

Reconstructing those events where the pions were emitted into the backward region by 
missing mass techniques in principle caused no problems. In several cases, however, we found 
events where three protons were detected as for a true dp —>■ pp7T~ps reaction, only there 
a third proton was produced in a chain of two consecutive np elastic scattering processes. 
From the first quasielastic scattering reaction one gets a scattered proton, a forward flying 
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Figure 2: Experimentally deduced momentum distributions for the spectator protons (solid 
histograms) and the reaction protons (dashed histograms) for two ranges in Q. 



spectator and a neutron. The scattered neutron in traversing one of the start detector 
elements hits another proton which reaches the detector whereas the slowed-down neutron 
remains unobserved. In simulating this process we found that by suitable selections in 
missing mass and angles these events could be eliminated. Thus an additional set of roughly 
0.6-10^ reconstructed ppir'ps events was obtained. 

As has been outlined in [12], the timing signals deduced from both ends of the Barrel scin- 
tillators not only yield information on the flight times, but also on the hit position of any 
track passing through the Barrel. Hence an important step in the detector calibration is 
the fixing of the absolute time offset which was carried out through a comparison with the 
results obtained for dp elastic scattering. This binary reaction with its unique kinematics 
and sizeable cross section was repeatedly measured in separate runs with an adjusted trigger 
condition. As a check of the reliability of the event reconstruction we show in Fig. 3 the 
deuteron angular distribution (given as histograms) in comparison with older data obtained 
at a somewhat higher beam momentum (solid dots) [H] . Instead of the deuteron beam mo- 
mentum of 1.85 GeV/c we quote a value of 0.92 GeV/c (T^^ = 376 MeV) which corresponds 
to the inverse reaction for a proton beam hitting a deuterium target at the same ^/s. The 
absence of data in the forward region is due to the fact that the corresponding protons were 
emitted towards angles > 60° which is out of the acceptance of the spectrometer. The over- 
all agreement is very good, the apparent mismatch in the peaking of the forward maximum 
results from the difference in beam momentum pp. To eliminate this dependency on pp we 
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Figure 3: Angular distributions for the elastic pd scattering reaction plotted as a function 
of the CM angle cos 9*^ (top) and momentum transfer t (bottom) in comparison with data 
from ref.^l4l- The dashed vertical lines denote the acceptance limits of our detector. 



plotted the same data (for cos 6^ > —0.35) as a function of the Mandelstam variable t 
(Fig. 3, lower frame) and found a very satisfying agreement. 



2.2 Monte Carlo simulation 

The analysis of our experimental data samples was accompanied by extensive Monte Carlo 
simulations. In order to allow each simulated quasi-free np pp7r~-event to have different 
initial kinematical parameters the program package was modified in a way as was described 
in detail in [12|, hence we will give only a short outline of the main ideas. Using the 
CERNLIB event generator GENBOD [15] one generates A^-body events for a given reaction 
specified by A^, type and mass of the particles involved and the total CM energy ^/s. The 
code returns momentum 4- vectors for each ejectile in the overall center-of-mass system and 
weight factors We based on the phase space density of the reaction. In the present case 
the basic reaction to be simulated is np ppn' . For each event randomly chosen values 
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for cos 9*, 0* and momentum |p*|/(MeV-c~^) were picked, the two former ones following 
uniform distributions, whereas the momentum was folded with the above mentioned Fermi 
distribution. We identify the three-component vector cos 9*,(f)* as well as the one 
pointing into the opposite direction with those of an np-pair within the deuteron in its 
CM system. Transformation into the laboratory system then allows one to deduce the 
corresponding vectors for spectator and projectile particle within a fast moving deuteron of 
momentum p^^ 1.85GeV/c. The fact that in the laboratory system the flight direction of 
the projectile neutron deviates by a small angle from that of the beam deuteron is accounted 
for by a suitably chosen rotation such that the neutron's flight direction serves as the actual 
beam direction. After having fixed event-by-event the momentum vector for the "beam 
neutron" it is straightforward to perform the simulation for np ppir' . 
By using approximately 1 million Monte Carlo events uniformly distributed across the avail- 
able phase-space we could determine the energy dependent acceptance of our detector and 
the reconstruction efficiency as a function of excess energy Q. The main limitations in ac- 
ceptance stemmed from the maximum in detector angle 6max = 60° and from the charged 
particles' energy loss in the various detector layers resulting in a low /3-threshold of j3 ~0.5 
for TT^ mesons and j3 ~ 0.35 for protons. In Fig. 4 we show the resulting acceptance curves 
for the relative proton momentum angle cos 9p (left) and the proton-proton invariant mass 
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Figure 4: Simulated acceptance curves as obtained for the relative proton momentum angle 
cos 9p (left) and the proton-proton invariant mass Mpp (right) for different values of excess 
energy Q. The dashed vertical lines in the plot of the right panel denote the kinematical 
limits. 
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Mpp (right) for selected values of Q. The acceptance goes to zero near |cos 9p\ = 1. This 
comes as a result of the way the relative proton momentum vector P* is constructed (see 
also Fig. 1). When the direction of P* approaches the beam direction, one of the protons in 
the CM system moves backwards, hence has minimum energy and drops below the detection 
threshold. Similar calculations have been performed for all other observables. 

3 Results and discussion 

Cross sections for various emission angles and invariant masses of each two-particle subsys- 
tem as well as two-dimensional Dalitz plots were extracted from the data. In order to derive 
absolute cross sections one must know the integrated luminosity. Defined a.s C = J rih ■ ritdt 
with nh{nt) denoting the number of beam and target particles, respectively, one finds the 
cross section from the relation 

a = n/(£-/-e), (2) 

where n is the number of observed events, / the deadtime correction factor and e gives the 
geometrical and reconstruction efficiency. This simple relation, however, has to be modified 
in case of a quasifree reaction. The Fermi motion of the neutrons within the deuteron will 
lead to a wide span in excess energy Q such that the number of beam particles initiating a pn 
reaction at a given Q will vary. In the present case of a close-to-threshold measurement, one 
furthermore will observe a strong variation in cr. In order to extract the energy dependence 
of the cross section and to compare it with the one of the free reaction it is necessary to 
unfold the effect of the Fermi motion from the data. By dividing the range in Q in small 
bins < Q >i such that per bin the variation in a is small and can be approximated by a 
constant o"j the number of produced events Ni is [16] 

N, = a,c[ |0(pb)M'p6, (3) 

J<Q>i 

where the integral is taken over all neutron beam momenta contributing to < Q >i and 
(p{pb) is the deuteron wave function as given by the PARIS potential [;17]. Here C is again 
the overall luminosity, its Q-dependence is accounted for by the integral. Correspondingly 
the number of observed events is given by 

n, = N,-f-e, (4) 
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The evaluation of the integral is performed by means of Monte Carlo simulations. Denoting 
the total number of generated Monte Carlo events by A^^"^*^ and the one generated for the 
bin < Q >j by N^'^^ the integral is given by the ratio 



Defining ni = ni/ei as the number of observed and acceptance corrected events the cross 
section is essentially given as cfj oc since £, / and N^^ are constants. This is 

demonstrated in Fig. 5 where we show the distribution of observed events in the top frame 
as the solid histogram; it extends from threshold up to 90 MeV. Also shown in this frame 
is the corresponding distribution as obtained for our Monte Carlo events. Calculated for a 
deuteron beam momentum of Pd = 1.85GeV/c it has its largest Q values also near 90 MeV, 
but on the low side starts with a sizeable yield already at threshold. Its maximum is shifted 
to lower Q values towards the peak of the deuteron wavefunction. When extracting the 
ratio of the experimentally deduced distribution and the Monte Carlo data, the histogram 
as shown in the bottom frame is obtained which (note the logarithmic scale) rises by more 
than two orders of magnitude. Also shown as a dashed curve are the total cross section data 
obtained with a free neutron beam at PSI [6] and parameterized as a S^'^ order polynomial 
in Q. When applying a suitably chosen normalization factor the present data are in good 
agreement with these absolute values with only some minor deviations at the lower and 
upper ends of the covered range. Henceforth this one normalization factor will be used in 
all of our further presentations and discussions of differential cross sections. 
No attempt was made to derive in an independent way absolute cross sections from the 
present experiment. The natural choice for the determination of the luminosity C would 
have been the pd elastic scattering reaction. Although it was quite successfully used for 
calibration procedures we did not consider it as being suited for finding the size of C Firstly 
the amount of available cross section data is still scarce. Apart from the already mentioned 
experiment by Booth et al. [14J at pp = 0.99GeV/c which corresponds to a bombarding 
energy of 425 MeV we only found one more set of published data by Alder et al. at 
comparable energies. These authors present data at proton bombarding energies of 316 and 
364 MeV covering far backward angles cos 9* < —0.6 and at 470 and 590 MeV at angles 
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Finally, by using eqs. [3-5] one finds for the cross section 
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Figure 5: Top: Simulated (dashed, MC) and measured (solid histogram, exp) distributions 
of pp7c~ events as a function of Q. Bottom: Ratio of the two distributions exp/MC (solid 
histogram) normalized to the cross section data measured at PSI 'i6] (dashed curve). 

cos 9* < 0. Due to this small range in cos 9* we consider these data unfit for a reliable 
interpolation. Recent data by Giilmez et al. |19| were taken at much higher energies of 
641 and 793 MeV, those by RohdjeB et al. [20] at energies up to 300 MeV. Secondly we 
found it difficult to estimate the error in cTpd when extracting the dp elastic events from the 
underlying background which was dominated by the much stronger pp quasielastic scattering 
events. Finally the uncertainties due to effects like shadowing and rescattering which tend 
to reduce the cross sections of any quasifree reaction by about 8% |21] and thus add to the 
size of the systematic error, would only then be of minor consequence when a comparison 
with another quasifree reaction is carried out. 

An order of magnitude estimate of the cross section could nevertheless be made. The inte- 
grated luminosity was found from the known target thickness (4 mm liquid hydrogen corre- 
sponding to 1.8 ■ 10^^/ cm^ target particles), the average beam intensity of 7-10^ deuterons/ s, 
and the total running time to be of order 40nb~^. Using eq. 2 with / ^ 0.5 and e ~ 0.25 
one finds a mean cross section near 60 /ib in good agreement with the PSI data [6]. 
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3.1 Angular distributions 

Acceptance corrected angular distributions are shown in Figs. 6 and 7 together with fits in 
terms of Legendre polynomials. The excess energy range 1.0-88.0 MeV was cut into six bins, 
namely [1.0-18.0], [18.0-34.0 ], [34.0-47.5], [47.5.0-61.0], [ 61.0-74.5], and [74.5.0-88.0] MeV, 
the indicated excess energies < Q > denote the center values of these bins. Error bars when 
given denote statistical errors only. It should be kept in mind that, although the cross 
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Figure 6: Angular distributions of the relative proton momentum for selected excess energies 
< Q> together with results of Legendre polynomial fits. 



section rises monotonically with Q, the observed counting rate does not. Due to the non- 
uniform Fermi distribution which governs the available excess energies, the highest rates 
are found near Q = 50 MeV and consequently one also observes the lowest statistical errors 
there. As already outlined above (see Fig. 1) the angular distributions of the relative proton 
momentum by construction are symmetric with respect to cos 9*p = 0. They are plotted as a 
function of cos^^^p and were fitted with even Legendre polynomials W{cos 6') oc 1 + ■ 
P2v{cos 6*) up to 1/ — 2. The extracted expansion coefficients are given in Table 1, up to 
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< Q > = 68 MeV the P4-term was neglected. In addition we give the numerical values for 
da /(El* in Table 2 where, as mentioned before, the absolute scale was adjusted to the PSI 
data [6]. 

Table 1: Expansion coefficients of the Legendre polynomial fits to the angular distributions 
of the relative proton momentum. 



< Q > /MeV 02 



10 


0.34±0.08 


26 


0.43±0.07 


40 


0.51±0.06 


54 


0.44±0.05 


68 


0.54±0.05 


82 


1.28±0.13 0.37±0.13 



Table 2: Differential cross sections in fib/sr for the angular distributions of the relative 
proton momentum at six excess energies Q. 





10 MeV 


26 MeV 


40 MeV 


54 MeV 


68 MeV 


82 MeV 


0.002 


0.124±0.022 


0.94±0.14 


2.04±0.16 


4.42±0.34 


7.96±0.61 


8.58±1.31 


0.014 


0.147±0.023 


0.87±0.13 


2.01±0.15 


4.50±0.36 


7.56±0.62 


8.78±1.43 


0.040 


0.142±0.024 


0.89±0.15 


2.07±0.16 


4.45±0.36 


7.69±0.64 


8.93±1.52 


0.078 


0.134±0.024 


0.92±0.16 


2.16±0.17 


4.67±0.38 


7.81±0.65 


9.29±1.48 


0.130 


0.150±0.024 


0.96±0.16 


2.26±0.20 


4.87±0.38 


8.30±0.71 


9.70±1.53 


0.194 


0.151±0.028 


0.98±0.16 


2.40±0.25 


5.10±0.46 


8.85±0.73 


11.07±1.57 


0.270 


0.154±0.033 


1.08±0.21 


2.60±0.33 


5.45±0.52 


9.65±0.92 


12.24±1.60 


0.360 


0.170±0.044 


1.18±0.34 


2.74±0.38 


6.04±0.54 


10.78±1.13 


15.21±2.13 


0.462 


0.120±0.047 


1.02±0.46 


3.20±0.52 


5.42±0.80 


11.33±1.31 


18.98±2.76 


0.578 


0.222±0.071 


0.90±0.61 


3.85±0.63 


6.15±1.11 


10.98±1.53 


21.94±3.12 



As one can see from inspection of Fig. 6 the scatter of the data points as well as the size of 
the error bars increases drastically for values cos^6'p > 0.4. This is the result of the very low 
acceptance observed in this angular region (see also the discussion in context with Fig. 4). 
Accordingly the minimisation was only performed on the first eight data points. In the 
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literature we found one measurement of these proton distributions, which was extracted 
from roughly 4000 bubble chamber frames [22]. At an average value of Q near 54MeV 
the authors report a value of 02 = 0.276 ± 0.032 which is to be compared to the present 
one of a2 = 0.44 ± 0.05. We believe the observed Aa deviation to be due to systematic 
errors in their method which were not included in the quoted error. As mentioned above 
only Ss, Sp, Ps and Pp partial waves should be present in the energy region covered in the 
present experiment. From our data, however, it can be seen that in the higher Q range 
contributions of the Ds wave are present as well. Near < Q > = 82MeV a P4(cos 6) term 
with a sizeable expansion coefficient 04 = 0.37 ± 0.13 had to be included in the fit. 
The pion angular distributions as deduced for the same intervals in Q are shown in Fig. 7. In 
general all are asymmetric and were fitted with sizeable ai and 02 coefficients (see Table 3). 
The one obtained at < Q > = 54 MeV is compared with data taken from ref. [22] (dotted 
line) and ref. [6j (dashed line) and good agreement is observed. The cross section as given 
in [22] exceeds the one measured at PSI by the factor 1.29. The authors of ref. |6j explain 
this discrepancy with a possible underestimation of the mean neutron energy in the older 
experiment. That measurement had been carried out over a broad neutron energy range 
and the mean energy had been deduced from a maximum likelihood fit. In the present 
comparison the data of [22] have been rescaled to match the PSI cross section. For the sake 
of completeness we additionally give in Table 4 the numerical values of the differential vr" 
cross sections. In passing we like to add that the corresponding ones given in ref. [6] (Table 
3 and Fig. 11) are not consistent with the absolute cross section data presented in their 
Table 4, but are too low by the factor 27r/10 due to an error in binning the data |23j . 

Table 3: Expansion coefficients of the Legendre polynomial fits to the it angular distributions. 
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Figure 7: { Color online) Angular distributions of the pions for selected excess energies <Q> 
together with results of Legendre polynomial fits. At <Q>=54 MeV, the pion angular 
distributions as found by Daum et al. 16] and Handler 122] are also shown as a dashed 
(blue) and a dotted (red) curve, respectively. 

3.2 Dalitz plots and invariant mass distributions 

Dalitz plots and Mpp invariant mass distributions of acceptance corrected and kinematically 
fitted ppTT~ events are presented in Figs. 8 and 9 together with statistical errors in the latter 
figure. In each case four Q bins 2 MeV wide were chosen, the center values are indicated 
in each frame. The kinematical limits of the Dalitz plots given by the solid lines were 
calculated for these values; due to the rapidly growing phase space some data extend over 
these border lines. Here the size of the squares is a measure of the count rate. Each plot 
is almost uniformly covered with the exception of the area in the upper left corner where 
strong FSI effects between the reaction protons were expected. Also some lowering in yield 
is observed in the opposite corner which we attribute to the asymmetries found in the pion 
angular distributions. No enhancements due to the A resonance are visible. 
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Table 4: Differential cross sections in /ib/sr of the pion angular distributions at six excess 
energies Q. 
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The da/dMpp invariant mass distributions shown in Fig. 9 were plotted on a linear Mpp scale, 
also given are the results obtained from our Monte Carlo simulation (solid and dashed lines). 
In all cases large deviations are observed between experimental and simulated data, as long 
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Figure 8: Experimentally deduced Dalitz plots for the quasifree reaction np — * ppn at four 
2MeV wide Q bins. The solid lines denote the kinematical limits. 

as purely phase space distributed events were considered (dashed hues). Incorporating FSI 
effects into our MC simulations by using the formalism of Watson [2l] and Migdal [25] which 
was later refined by Morton [26] the distributions given by the solid lines were found. We 
calculated additional weight factors Wfsi given in a simplified form as 



2 

1-1 



(7) 



where T^*^ denotes the pp center of mass kinetic energy Tg^ = {Mpp - 2mp) and the 
Coulomb penetration factor 



with 7p = "'^pp''^ , Here a is the fine structure constant, Ppp = y'2/ippT^^^ and fipp is the 
reduced mass of the pp-sjstem. The strength factor fpp is a measure of the contributing 
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Figure 9: (Color online) Proton-proton invariant mass distributions obtained at four Q bins 
together with data from our Monte Carlo simulation. The dashed (blue) lines denote the 
results as found for phase space distributed events, the solid (red) curves give the ones where 
FSI effects with standard values for scattering length and effective range have additionally 
been included (see text). 



Ss and Sp partial waves and is adjusted for each Q interval. From literature we took the 
standard values ao=-7.83 fm and ro=2.8 fm [27J as input parameters for the scattering 
length and effective range, respectively, for the two protons in the ^So state. The agreement 
for the two lowest Q bins, where the relative weight of the "diproton" IS'o-state is high, 
is very good. In case of the two higher bins this simple ansatz, however, only succeeds in 
reproducing the rise at Mpp = 2mp. 
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3.3 The isoscalar cross section 

Using equation (1) the isoscalar cross section aoi can be obtained from the measured cross 
sections for the pp ppn^ and the np ppTT~ reactions. The isospin 1 = partial waves 
of type Sp (^5*1 ^^SqPi and ^Di -^^SqPi) which are forbidden in the pp — > ppn'^ reaction 
due to the Pauli principle are generally believed [2,6,10,11] to dominate cxoi in the threshold 
region. In the np reaction they interfere with the isospin 1 = 1 wave ^Pq -^^SqSq and thus 
are responsible for the strong asymmetries in the vr" angular distributions at very low 
excess energies. In the literature one finds a vast amount of data for both reactions and, 
as has been shown [6], the extracted cxoi shows the expected t]^ dependence at least for 
rj > 0.5 {Q > 17MeV). We contend that the deviations quoted for smaller rj values are 
the result of wrong cross section data. As can be seen from (1), in order to obtain a finite 
ctqi the np cross section must at least be half as large as the one for the pp reaction. For 
T] = 0.34 (T„p is given as 1.43 fib [6]. Recently our group at COSY reported new data for the 
pp ppiT^ reaction which exceeded the published ones from lUCF [28] and CELSIUS [29] 
by roughly 50%, the discrepancy could be shown to originate from an underestimation of 
the pp final-state interaction [30] • We found a cross section cipp = 3.72±0.3/i6 at ?7 = 0.35 
which is 2.6 times larger than cr„p and as such would leave no room for ctqi- This, however, 
is in contradiction to the asymmetries observed for the tt~ angular distributions, which are 
only possible with interfering Ss and Sp partial waves. In addition to assuming a wrong 
cross section measurement one should also consider a wrong beam energy determination. 
The uncertainty in neutron energy at the NA2 beam facility at PSI is given as 3 MeV (a) 
for En = 287 MeV [6|. An error in quoted beam energy of this size could possibly explain the 
observed deviations in the very close to threshold region where the cross section of almost 
any reaction rises dramatically. 

3.4 Summary 

With a deuteron beam at 1.85GeV/c impinging on a liquid hydrogen target the quasi- 
free np pp7T~ reaction was studied for excess energies Q up to 90 MeV. The data were 
analyzed in the framework of the spectator model where the proton is assumed to be an 
unaffected spectator staying on-shell throughout the reaction. Tagging the spectator proton 
in the forward scintillator hodoscope of our COSY-TOF spectrometer allowed to determine 
such parameters as effective mass and momentum of the off-shell neutron at the time of the 
reaction. We have measured angular distributions and invariant mass distributions of the 
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reaction products and have set up Dalitz plots for several Q bins distributed evenly over the 
whole excess energy range. The data were compared to results derived from Monte Carlo 
simulations and to data taken from the literature. In general good agreement was found, the 
large asymmetries observed previously for the vr" angular distributions could be confirmed. 
Final-state interaction effects between the reaction protons were found even at the highest 
excess energies. The da/dMpp invariant mass distributions at 25 and 40MeV which are 
governed by the "diproton" ^S'o-state could be reproduced by the Monte Carlo simulations 
when incorporating FSI effects in the formalism of refs. [23-25] with standard values for 
scattering length and effective range. Sizeable d-wave contributions were observed in the 
angular distribution of the relative proton momentum at Q = 82 MeV. In view of new cross 
section data of the pp —>■ ppn^ reaction, reported deviations of the isoscalar cross section ctqi 
from an rj^ dependence [B] were explained as stemming most probably from small errors in 
beam energy. 
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